A simple microscopic description of quantum Hall transition without Landau levels 
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By restricting the motion of high-mobility 2D electron gas to a network of channels with smooth 
confinement, we were able to trace, both classically and quantum-mechanically, the interplay of 
backscattering, and of the bending action of a weak magnetic field. Backscattering limits the mobil- 
ity, while bending initiates quantization of the Hall conductivity. We demonstrate that, in restricted 
geometry, electron motion reduces to two Chalker-Coddington networks, with opposite directions of 
propagation along the links, which are weakly coupled by disorder. Interplay of backscattering and 
bending results in the quantum Hall transition in a non-quantizing magnetic field, which decreases 
with increasing mobility. This is in accord with scenario of floating up delocalized states. 

PACS numbers: 72.15.Rn; 73.20.Fz; 73.43.-f 



Introduction. Quantization of the Hall conductivity of 
a disordered 2D electron gas, axy — n, (in the units of 

/h) together with vanishing diagonal conductivity, axx, 
reflect the fact that in a perpendicular magnetic field 
delocalized states always constitute a discrete set [l;]. 

In a strong magnetic field, uJcT ^ 1, where tUc is a 
cyclotron frequency and r is the scattering time, energy 
positions, of the delocalized states coincide with the 
centers of well-resolved Landau bands. Such a strong- 
field limit was the focus of theoretical studies of delocal- 
ization in a magnetic field. Most appealing qualitative 
picture Q assumes a smooth disorder when the eigen- 
states are well-defined Larmour circles drifting along 
equipotential lines. Then delocalization corresponds to 
the classical percolation threshold; localized states above 
En are closed drift trajectories executed, e.g., clockwise, 
while the states below £"„ are closed drift trajectories 
executed counter-clockwise, see Fig. [1] 

An alternative approach ^3] to delocalization is based 
on renormalization-group equations, describing the evo- 
lution of (Jxx, cFxy upon increasing the sample size, L, 
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where _D is a dimensionless constant. First term of 
Eq. ([T]) originates from interference of electron multiple- 
scattering paths: two paths corresponding to the same 
scatterers but different sequences of scattering events in- 
terfere even in the presence of Aharonov-Bohm phases. 
Second term reflects the orbital action of magnetic field: 
by curving electron trajectories it tends to destroy the in- 
terference. When the "phase" and "orbital" terms com- 
pensate each other, delocalization transition takes place. 

Field-theoretical approach H yields a highly nontrivial 
prediction first pointed out by Khmelnitskii [4l| . Namely, 
solving Eqs. ((T|), ^ together with classical initial con- 
dition axy{i^c) = c!QbJcT{\ -f oj^r^)"-'^, where (Tq oc En is 
the dimensionless conductance at ojc = 0, yields En = 



huc {n + ^) [l + (uJct)^^] . As shown in Fig.[T]for n — 0, 
the high-field part, oJcT ^ 1, oi Eq follows the center of 
the lowest Landau level, while the low-field part "floats 
up" as ujcT — > 0. This prediction is essential component 
of the global phase diagram ! 5| . 

Qualitative classical picture 2] applies to the high-field 
part and illustrates the restructuring, see Fig. [TJ of the 
motion of guiding center, which accompanies the crossing 
of Eq by the Fermi level, upon increasing magnetic field 
(1 — > transition into the quantum Hall insulator tel). 
While there is certain experimental evidence 0, H, H, 



3, El, E, El that floating up of En 



indeed takes 

place, tight-binding numerical studies [14l.ll5l. 16l| are less 
conclusive. There is a fundamental reason [l7| for this 
lack of conclusiveness. Indeed, significant fioating up oc- 
curs for large ctq > (wct)~^ ^ 1. For such ctq, upon mov- 
ing along the dashed horizontal line in Fig. [1] the local- 
ization length undergoes a steep growth from orthogonal 
value exp(7rtTo) to the unitary value ^„ exp(7r^(TQ) 
already in very small fields uJcT ~ aQ^e~'^^'^° , way before 
the expected transition at point L. Such large is a ma- 
jor obstacle for numerics to capture the descending region 
of the line of transitions in Fig. [TJ Accessing this region 
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FIG. 1: (Color online) Energy position of the lowest delocal- 
ized state, Eo, vs. magnetic field, lUc, as predicted in Ref. Q. 
For a given E, delocalization transition occurs at points L 
and H. Change of character of electron motion near H is 
illustrated schematically. Corresponding change near L is il- 
lustrated in Fig. [5](a «^ e). 
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requires to construct an efficient minimal description of 
a weak-field transition, as transparent as the picture 
sketched near the point H in Fig.[TJ This goal is achieved 
in the present paper. The key step of our construction 
is separation of the spatial regions with disorder-induced 
scattering and field-induced bending. 
Restricted electron motion. (i) In contrast to uni- 
directional motion in strong fields, we allow counter- 
propagating paths in the regions where orbital action of 
magnetic field is negligible. We achieve this by restrict- 
ing electrons to narrow point contacts, see Fig. [D At the 
same time, we assume that the phase action of magnetic 
field is well-developed in each point contact, i.e., the area 
of the contact is threaded by many flux quanta. Presence 
of disorder is incorporated by allowing mutual backscat- 
tering of two counterpropagating waves. We quantify the 
strength of backscattering with probability, p, so that the 
scattering matrix of the contact has the form 

with amplitudes Zi, Zi, Fig. [21 having random phases. 
{ii) The orbital action of magnetic field takes place in 
the junctions between the point contacts, Fig. [2] To 
simplify the description of the junction, we assume that 
an electron incident, say, from the left, after several 
bounces [3| off the walls exits either "up" or "down", 
i.e., both forward and backward scattering channels are 
suppressed. This assumption allows us to quantify the 
bending strength of the junction by a single parame- 
ter, q, the deflection probability to the right, Fig. [21 
Then the deflection probability to the left is {1 — q). 
The Hall resistivity of the junction [13] is then given by 
Rh = ieyh)[q - (1 - + (1 _ q)% so that q = 1/2 

corresponds to a zero field. Expression for the scattering 
matrix of the junction is the following 
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With scattering matrices Eqs. (HI) defined, the prob- 
lem of electron localization by disorder in a magnetic field 
reduces to the effective network model, which can be 
studied by transfer-matrix method, similar to Chalker- 
Coddington (CC) model which describes delocal- 
ization transition at the point H in Fig. [1] As in Ref. 
[igl . "unitary" disorder is incorporated via random phases 
of the link amplitudes, Zi. Delocahzation transitions in 
the network define a line on the (j), q) plane. Important 
is that this line can be converted into the dependence 
i?o(wc). Indeed, parameter q refiects the strength of mag- 
netic field, so that (5 — 9) oc cjc, while the backscatter- 
ing probability, p, decreases monotonously with increas- 
ing energy. Thus, the floating scenario is equivalent to 




FIG. 2; (Color online) Left: Restricted electron motion over 
point contacts and bend-junctions is illustrated; An,m are the 
centers of forbidden regions. Right: Scattering matrices of 
the junction and of the point contact. 



the statement that p(g)-line approaches p = as g ap- 
proaches i. Below we argue that the form of p(g)-line of 
delocahzation transitions is the one shown in Fig. [S^a) 
(region q < i), so that it indeed yields the dependence, 
p{q), corresponding to the floating of £^o(^c)- 

In the CC model the transmission of the nodes with 
"hei ght" , Ei, and "width", F, is given by the Fermi func- 
tion ;20] T{E,E,) = {1 + exp [{E, - E)/T]}~^. Qualita- 
tive strong-field picture of the transition Q emerges when 
the spread, W, of heights, Ei, is ^ F. Then the quan- 
tum interference can be neglected up to large distances, 
Ti-ci = {W/T)'^^'^ 1, determined by the classical perco- 
lation. For smaller distances, one can replace T{E,Ei) 
by a step-function, Q{Ei — E). Adopting the same ap- 
proach, we assume that (i) in Fig. [21 full transmission 
takes place in p percent of point contacts, and full re- 
fiection in the rest (1 — p) percent; (ii) a given junction 
deflects only to the left in q^ percent of cases, only to 
the right in (1 — g)^ percent of cases; in the remaining 
2q{l — q) percent the deflection takes place both to the 
left and to the right depending on incoming channel. 




FIG. 3: (Color online) (a) Classical metallic phase, dashed 
region (e) in the {p, (jr)-plane separates different insulating 
phases of the model Fig. [2] (b) Line of classical delocahza- 
tion transitions in the (p, g)-plane. Asymptotes Eqs. (??), 
([B]), are shown with green and black lines, respectively; red 
curve interpolates between the two regimes. Crosses show the 
positions of quantum delocahzation transition inferred from 
simulations for five values of "energies", p. Quantized axy is 
1 in phase (a), —1 in phase (c), and in (b), (d), and (e). 
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Phase diagram. The key observation that allows to estab- 
lish the phase diagram Fig. [3]is that the classical electron 
motion over the lattice of point contacts and junctions 
can be reduced to a single problem of joint bond percola- 
tion over "p" and "q" -bonds. To substantiate this state- 
ment, we focus on the grey squares in Fig. [21 which are 
"forbidden" regions for electrons, and notice that elec- 
tron scattering processes both in point contacts and junc- 
tions effectively establish bonds between these regions. 
More specifically, if electron is backscattered in a point 



contact, we consider that the centers An m and A„ 



-lOf 



the squares, adjacent to this contact, are connected by a 
bond, see Fig.|4]. Further, if electron is bent-scattered by 
a junction, say, in the direction left — > down, we identify 
this process with establishing a bond between the centers 
of the squares An^m and 74„_i,„i_i. The above identifi- 
cation reduces the classical motion through the network 
with sites in the centers of squares, see Fig. ID 
Structure of phases. We start from the region of strong 
reflection, p > ^. Counterpropagating channels in point 
contacts, Fig. [2J are essentially "short-circuited" . Hence, 
no delocalization occurs upon increasing magnetic field, 
(^— g). Localized states are illustrated schematically 
in Fig. [5] in the limit {1 — p) ^ 1. It is apparent that 
crossover from clockwise rotation (b) to counterclockwise 
rotation (d) upon passing the zero-field line q = takes 
place without delocalization. This absence of delocaliza- 
tion is consistent with low-energy part of Fig. [1] because 
P > \ corresponds to small do- Along the line P — \ the 
p -bonds alone constitute a critical network. It is seen 
from Fig. [5] that, as p is reduced below i, q -bonds sus- 
tain percolation, manifesting that metallic behavior for 
(To > 1 persists up to strong magnetic fields, oJcT ^ 1. 

We now turn to the most interesting region of small p 
(high energies in Fig. [T]). In this domain the overall con- 
nectivity of the network is dominated by the q -bonds. 
Moreover, at p = the light-blue and dark-blue subnet- 
works are completely decoupled. Small finite p becomes 
essential in the vicinity oi q~ i , when both subnetworks 
of q -bonds are critical. Now even weak non-zero p, by 
coupling the subnetworks, results in opening of a metallic 
region (e) in the phase diagram. This disorder-induced 
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FIG. 4: (Color online) (a) In the classical limit, electron mo- 
tion reduces to the bond percolation on the lattice consisting 
of p and q -bonds; (6) and (c) are graphic representation of 
Eqs. (??) and (O, respectively. 



coupling of critical subnetworks is the fundamental un- 
derlying mechanism for the restructuring of states near 
the point L, Fig. [TJ In effect, transformation (e) — > 
(a) with increasing field, (5 ~ <?), is the counterpart of 
transformation near the point H in quantizing ujc- 

Fig. [5] (a) also illustrates that at small p both subnet- 
works are chiral. Transformation into the phase (e) upon 
decreasing magnetic field is accompanied by the change 
of the Hall conductivity from quantized to finite value 
smaller than 1 . Full suppression of axy in the region (e) 
occurs only when interference drives this region into the 
Anderson insulator, so that the difference between (e) 
and strongly localized phases (b) and (d) vanishes. 

Calculation results for the phase boundaries are shown 
in Fig. ^h) only for g < ^ due to g ^ (1 duality [2H . 
The end-points {p,q) = (O, ^) and {p,q) — (^,0) are 
conventional bond percolation thresholds for percolation 
over regions A^.n, connected by q -bonds and by p - 
bonds, respectively, see Fig. [H Next, consider small p ^ 
1 and suppose that q is slightly smaller than i. Then, 
instead of missing q -bond, a pair of one horizontal and 
one vertical p -bonds can provide the connection (there 
are two such variants, see Fig. [4^)). Resulting shift of the 
threshold position is determined by the condition 

/(P, q) = q' = q + 2(1- q)p\l + q) ^ 1/2, (5) 

where q' is the probability that effective q -bond connects. 

When q is small, the role of q -bonds is to promote 
percolation over p -bonds. Unlike the previous case, the 
shift, (i — p), of the threshold is linear in q. This is 
because in certain rare cases one q -bond takes on the 
role of a missing p -bond. Assume that a vertical p 
-bond is absent while one of the neighboring horizon- 
tal p -bonds is present. Then the connection between 
the ends of the missing bond can be established in two 
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FIG. 5: (Color online) Illustration of the structure of five 
different phases in Fig. [3] Electron moves predominantly 
along solid directed lines. In phases a and c the motion is 
either within dark or within light subnetwork. In phase e 
coupling of subnetworks by p -bonds opens a metallic phase. 
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steps: first via this horizontal p -bond and then via a 
q -bond, Fig. [4] Quantitatively, the boundary p{q) at 
small q can be obtained from the real-space renormaliza- 
tion group procedure . In Ref . [l^ the probability that 
the supcrbond, illustrated in Fig.[3J: connects, is given by 
Mp) =p^ + 5p'^(l-p) + 8p3(l-p)2 + 2p2(i-p)3^ where 

the last three terms correspond to realizations when su- 
perbond connects with one, two, and three original bonds 
removed; fo{p) — \ yields the exact threshold p = i. Us- 
ing the fact that q -bonds can restore the connectivity, 
and selecting suitable realizations out of all 2^ possible 
states of q -bonds amounts to the following modification 
of the probability, p' = f(p, q), that superbond connects 



length, £,M{p,q), was inferred from the smallest positive 
exponent: (,m{p,q) = ^~m/2- Simulation confidently con- 
firm that for finite "energies" , p > 0, quantum system 
is insulating at g = 0.5, while the state with p = and 
q = 0.5 is extended. As seen in Fig. [31d, the discrepancy 
between classical and quantum treatments is small. 

Note in conclusion that among various network models 
studied 25 1, the closest to ours is the model [1^. Unlike 
Ref. [Hour Eq. (|4]) describes scattering, say, to the right, 
with the same probability. 

Acknowledgments. This work was supported by the BSF 
grant No. 2006201. 



f{p,q) = fo{p) 

+p(l-p)4[4g- 
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-Sq^ + q']+{l-pf[2q'-q']. 



Upon equating f{p,q) to i, Eq. ^ yields the boundary 
of percolation transition at small q. 
Localization length. The boundary a-e of "classical" 
phase diagram Fig. [3l is characterized by the critical 
exponent v = ^. Simulations, see below, indicate that, 
with quantum interference, the phase diagram Fig. re- 
mains unchanged, while metallic phase in Fig. [5^ turns 
into the Anderson insulator with "unitary" localization 
length, ln^„ = tt^ctq, where the Drude conductance is 
related to p as p — ca^^ , with c ^ 1. Interference also 
modifies the divergence of the localization length 

Hp, q) - \f{p, q) - l/2r = |3pV2 - (1/2 - q) r^ (7) 

in the neighboring insulating region. Fig. [5^, by changing 
1/ from I' = |, to quantum « ^. Starting from small 
field, — q) ^ P^, S, retains its value ^„ oc exp(7r^c^/p^) 
up to the narrow vicinity of the boundary, ~ ?) = §P^, 
when it crosses over to diverging Eq. ([7]). The width of 
the "quantum" region can be estimated as 
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This region rapidly narrows in course of floating up. 
Quantum treatment of the network. Numerical simula- 
tions, employing matrices Eqs. ([3]), ^ for nodes and 
incorporating random phases into the link amplitudes, 
are required to verify the above predictions based on the 
classical picture. They are also supposed to verify that 
domain (e) in Fig. [3] is, quantum-mechanically, insulat- 
ing. The results of transfer-matrix analysis of the two- 
channel [2^ network Fig. [2] are shown in Fig. [Sb for 
five values of "energy", p. Usual simulation procedure 
2J] was employed: upon constructing a transfer matrix 
of a slice (M nodes in transverse direction with periodic 
boundary conditions) the net transfer matrix of a system 
of length, N (typical N = 240000), was obtained and 
diagonalized, yielding the Lyapunov exponents, \i{p,q), 
related to the eigenvalues as exp(AiiV). Localization 
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